Magnetic susceptibility in quasi one-dimensional Ba2V309: chain segmentation versus 

the staggered field effect 
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A pronounced Curie-like upturn of the magnetic susceptibility x{T) of the quasi one-dimensional 
spin chain compound Ba2V309 has been found recently Frequently this is taken as a signature 
for a staggered field mechanism due to the presence of g-factor anisotropy and Dzyaloshinskii- 
Moriya interaction. We calculate this contribution within a realistic structure of vanadium 3d- and 
oxygen 2p-orbitals and conclude that this mechanism is far too small to explain experimental results. 
We propose that the Curie term is rather due to a segmentation of spin chains caused by broken 
magnetic bonds which leads to uncompensated S=| spins of segments with odd numbers of spins. 
Using the finite-temperature Lanczos method we calculate their effective moment and show that 
~ 1 % of broken magnetic bonds is sufficient to reproduce the anomalous low-T behavior of xiT) in 
BaaVsOg. 



PACS numbers: 75.10.Jm, 75.40.Cx, 75.50.Ee 



I. INTRODUCTION 

Recently a quasi one-dimensional magnetic behavior of 
Ba2V30g was clearly revealed by means of the magnetic 
susceptibihty x(T) and the specific heat Cp(T) measure- 
ments in polycrystalline samples 0. The data are com- 
patible with the spin S = 1/2 antiferromagnetic (AF) 
Heisenberg chain model with a nearest neighbor exchange 
coupling J = 94 K. In addition an anomalous Curie-like 
upturn of xiT) was found below 20 K and claimed to 
be of intrinsic nature because the effect of paramagnetic 
impurities was ruled out by the analysis of experimental 
data. This low-T behavior of x(T) was tentatively at- 
tributed in U to the staggered field effect induced by the 
applied magnetic field. In a quasi one-dimensional spin 
S = 1/2 chain a low-T upturn of x(T) is expected if the 
Dzyaloshinskii-Moriya (DM) interaction and/or a stag- 
gered (7-factor anisotropy are present in the system 
Among the 3cZ-systems, Cu benzoate 5|| and pyrimidine 
Cu dinitrate are the most known examples. How- 
ever the anomalous low-T part of xiT) in Cu benzoate is 
much larger than in the theory 3] and the physical rea- 
son for this discrepancy is not yet clear. For better un- 
derstanding, a comparative analysis of different sources 
contributing to the anomalous low-T behavior of x(T) 
is necessary. In the present paper such an analysis is 
developed for Ba2V309 compound. 

The low symmetry of Ba2V309 allows both a stag- 
gered (7-factor anisotropy and a DM interaction to con- 
tribute to the staggered field. We suggest a simple model 
of vanadium 3d- and oxygen 2j)-orbital structure within 
magnetic chains and calculate both contributions. Alter- 
natively, we consider a segmented spin chain model where 
uncompensated spin moments lead to the low- T Curie 
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term and compare the results of both models to the ex- 
perimental xi^)- A segmented AF spin-1/2 chain model 
was applied earlier 0, Q to explain the magnetic suscep- 
tibility measurements in quasi one-dimensional cuprates. 
In these studies, analytic results |2J quantum Monte 
Carlo simulations j3| were used, assuming a particular 
random distribution of broken magnetic bonds (defects). 
Here we extend the segmented-chain model and com- 
pare different distributions of defects or chain segment 
lengths. To determine the thermodynamic properties 
of the chains under consideration, we apply the finite- 
temperature Lanczos method 



II. ORBITAL STRUCTURE AND THE 
3-TENSOR ANISOTROPY 

According to 0] VOe octahedra in Ba2V309 form 
edge-sharing chains. The Og octahedron with an aver- 
age V-0 distance ~ 2A is only slightly distorted, and 
V^"*" ion ((i^-state) is displaced from the center by about 
0.2 A towards one of the oxygens shared by two neigh- 
boring octahedra. The direction of V-ion off-center dis- 
placement distinguishes clearly a local z-axis in each VOg 
octahedron. The energy splitting of vanadium c?-orbitals 
was suggested Q to be such that c?2;j,-orbital located in 
the plane transverse to the local z-axis is the singly oc- 
cupied ground state orbital (Figure Results of ab 
initio calculations of a small cluster confirm this picture. 
Within a chain the local z-axis varies in a zig-zag manner 
thus forming alternating short and long V-0 bonds. The 
resulting c?-orbital arrangement along a chain determines 
the symmetry and strength of a superexchange (SE) cou- 
pling between vanadium spins. Based on standard SE 
theory we calculate both the isotropic exchange constant 
Jij and the DM vector for neighboring V-ions. 

First we calculate the g tensor components of the V''+ 
ion in the VOg octahedron whose size is scaled to the one 
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FIG. 1: Two edge-sharing VOe octahedra forming a segment 
of a magnetic chain in Ba2V309. The local coordinate sys- 
tems attached to the A- and B-octahedra and the directions 
of the off-center displacements of the VA,B-ions along the lo- 
cal z-axes are shown. The ground-state d^y and d^y orbitals 
of the VA,B-ions and the relevant p^i-orbital of a common- 
edge oxygen ion responsible for the Va-Vb superexchange 
are shown. The dxy orbitals in neighboring octahedra lie in 
orthogonal planes. The Dzyaloshinskii-Moriya vector Dab is 
parallel to a-a^ds and alternates in sign on successive bonds. 



in Ba2V3 09. Strong variation of V-0 covalent bonding 
due to a V-ion displacement from the central position in 
the nearly cubic Og cage splits t2g and eg vanadium d- 
orbitals. We neglect the extra distortion of the crystal 
field on V-ions caused by the side ionic groups attached 
to the chains. Then the diagonal g tensor is given by 
5|l = 2(1 - Azz),9± = 2(1 - A±), where = A^^^ = 
Ayy, and A^^ 
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^E,n^omL^.H\ViEm - Eo). Here, 
0.03 eV is the constant of spin-orbit coupling for 



V^^-ions is the orbital angular momentum and 

|0), |m) are ground state and an excited d-orbital state 
respectively with an excitation energy Em — Eq ^ X. 

For a VOe octahedron with one unpaired electron on 
vanadium d-shell wc performed ab initio quantum chem- 
ical calculations of the ground state Eo and excited state 
Em cluster energies as function of the ion displace- 
ment 5 along a local z-axis. The model cluster [VOg]*" 
was embedded in an octahedral point charge environ- 
ment and the program package molpro was used. First 
the restricted Hartree-Fock molecular orbital (MO) wave 
functions were determined to define an active space of 
MOs including a variety of valence and virtual orbitals 
(altogether 23 orbitals arising from strongly hybridized 
Y{3d) and 0{2p) states). Then the active space was 
used to take into account effects of electron correlations 
by performing multireference configurational interaction 
(MRCI) calculations of the cluster energies when the 
unpaired electron was placed successively in different d- 
orbitals. The results of these post Hartree-Fock calcu- 
lations (Figure ^ show that t2g-orbitals are split into 
the ground state orbital singlet \dxy) = |0) and degen- 
erate orbital doublet \dxz) = \dyz) = |2). At the 
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FIG. 2: Excitation energies Em — Eq of the Vanadium 3d 
states |m) in VOe as functions of the off-center displacement 
S of the V*+ -ion along the local octahedral z-axis. 



physical displacement S =0.2 A, the energy difference 
between them is Ei,2 — Eq ~ 0.23 eV = Ai, while for 
\dx2_y2 EE |3) one has E^ - Eq ~ 1.26 eV = A3. Us- 
ing (1|L^|0) = -(2|Lj,|0) = i and (3|L^|0) = -2i one 
obtains the following result for the (7-factor anisotropy: 
(311 - 5i) ~ 2A(1/Ai - 4/A3) ~ 0.07. 

Next we define a global reference system with unit vec- 
tors X II a, y II b and z || a x b , where a and b are or- 
thogonal crystallographic axes in Ba2V30g. In this sys- 
tem the g tensor is decomposed into a uniform and a 
staggered part, gi = g^ + (— l)^gs, where £ denotes the 
V-sites along a chain. The only non-zero components of 
gs are {gs)yz = {gs)zy = 9s = [g\\ - 5±) sinacosa, where 
a is a tilt angle. In Ba2V309 the local z-axes are tilted 
from the chain direction (b-axis) and a ~ 45° alternates 
as t runs along the chain. Finally, our calculation leads 
to a staggered component — 0.035 for Ba2V309 . 



III. ISOTROPIC EXCHANGE AND 
DZYALOSHINSKII-MORIYA INTERACTION 

Considering now the intrachain magnetic coupling we 
note that for each (AB)-pair of neighboring V-ions the 
shortest SE path is via the out-of-plane pj;-orbital which 
belongs to one bridging oxygen as shown in Figure ^ 
To estimate the electronic hopping parameter tdTx of the 
TT-bonding between each of occupied d^^^j^-orbitals and 
bridging oxygen p^j-orbital we use Harrison's prescrip- 
tion [l3: idTT = ripdirh^r^J'^ /{moJ/'^), where ripd^r = 1-36, 
f? jm, = 7.62 eVA and Vd = 0.98 A for vanadium. For the 
average V-0 distance a ~ 2 A one obtains id-K — 1 eV. 
Then from SE theory the isotropic exchange constant is 
obtained as 
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where Ud and 17^ are the on-site Coulomb integrals on 
V and O, respectively; A^^^ — id — ^tt is the energy dif- 
ference between the vanadium ground state d- and oxy- 
gen pj;-orbitals. We have A^^ > 0, since Ba2V309 is 
a Mott-Hubbard insulator. We use the following pa- 
rameter values (in eV): Ud — 3 ... 4, Utt — 2 ... 6, and 
AdTT — 5, characteristic of vanadium oxides with tetrava- 
lent V-ions. This leads to an exchange J in the range 
10 meV < J < 20 meV in comparison to the experimen- 
tal value J « 10 meV. We conclude that the use of the 
shortest SE path via the common edge 0-ion provides 
a satisfactory account of the isotropic AF exchange in 
BaaVsOg. 

For a pair (AB) of V-ions the DM vector can be 
written [ll| as Dab = — i(AA — Ab) with A^/b = 
2AEm5^o('^A/B|LA/B|0A/B)(^m " Eq)-^ j'pl- Here, 

= J(r7iA,0B;0A,0B) and J^'' = J(Oa, tob; Oa, Ob) 
are SE constants in similar intermediate configura- 
tions where the unpaired electron on the Va- or Vb- 
ion, respectively, is raised by the spin-orbit interac- 
tion to mth excited c?-state. The symmetry requires 

J(toa,0b;0a,0b) = J(0a,tob;0a,0b) = Jab^- In this 
notation, the isotropic constant calculated above reads 

J = 2 J(0a, Ob; Oa, Ob) = 2 J^g. For the orbital geometry 
in Figure ^ inspection shows that the only non-zero ex- 
cited state exchange Jab' m = 1. This corresponds 
to the intermediate electronic configuration with occu- 
pied \dxz) state. The matrix elements of La and Lb are 
calculated by referring to a common coordinate system, 
which yields (1a|La|0a) - (1b|Lb|0b) = -2ieAB, where 
the unit vector gab is defined as bab = za x zb- To com- 
pare J^g and Jab, we note that the hopping parameter 
t[}^ between each of the excited vanadium d^i^ -ovhitals 
and the bridging oxygen pa;-orbital is related to tdTi- 
hopping as t^^ ~ td-K sin/3, where t^Tr enters into the defi- 
nition of jf^ and sin/3 ~ (S/a) ~ O.l. This immediately 
leads us to the simple result, J^g « J^gSin/3, hence 
Dab/Jab = d-eAB, where d = (2A/Ai) sin /3 « 0.025. 

Dab is staggered along the chain direction, i.e. 
T>£/+i = (— 1)^D due to the property of e^^^+i. The geo- 
metric factor sin /3 is a measure of the asymmetry caused 
by off-center displacements of V-ions. If they are ne- 
glected the inversion center at the midpoint of (AB)-pair 
is restored and Dab = 0. 



IV. STAGGERED-FIELD MODEL FOR THE 
SUSCEPTIBILITY UPTURN IN BA2V3O9 

Taking both effects described above into account, a 
magnetic field H applied to a chain induces a staggered 
field h, which can approximately be written as 

h « ^D X H -f g,H. (2) 



For T <C J/kB, the staggered field leads to a contribution 
to the magnetic susceptibility of the form 3] : 

which scales with the factor c ~ h/H-^. This should be 
compared with the experimental low-T (2 K < T < 20 K) 
behavior of x in Ba2V30Q which is well described by a 
Curie law, xlt = Clt/tQ. 

The staggered susceptibility XsiT) varies approxi- 
mately like Cs/T. For 2K < r < 20K, we take into 
account the slowly varying logarithmic correction by re- 
placing ln^/^( J/A:bT) with its average value « 1.6. Then, 
we obtain the Curie constant due to the staggered field as 

a 0^ 1.33c2(7VAM0Ai|/3A:B) = (AfAMo/3fcB)(Mf )', where 
the second relation defines the effective staggered mag- 
netic moment = Vl.33c/LtB. 

In general the factor is an angle dependent func- 
tion of the magnetic field H direction: c^{9,(p) — 
J2fj,{9hfj,/dHa) , a = a{9,(f>). In a polycrystal used 
in 0, the angular average = (c^(0, 0)) is measured. 
With D and given, one obtains = (2/3) [d^ + g'^], 
yielding Cg — 2.5 • 10"'^ cm^K/mol, and correspondingly 
an effective moment fi^^ ~ 4 • 10~^/LtB. 

In comparison the experimental values reported in , 
namely Clt — 6.0-10"^ cm^K/mol and /i^x — 0.2/iB, are 
much higher. Therefore we conclude that the staggered- 
field effect alone is not sufficient to explain the low- 
temperature behavior of the magnetic susceptibility in 
Ba2V3 09, and a different mechanism must be present. 

V. SEGMENTED-CHAIN MODEL FOR THE 
ANOMALOUS SUSCEPTIBILITY 

The previous magnetic model implies a long-range or- 
der of short V~0 bonds in each structural V02-chain. 
From Figure n we infer that several quasi-degenerate en- 
ergy minima for the short V-0 bonds orientation can 
exist. This expected structural degeneracy is described 
in an easy way by recalling that a local z-axis in each VOg 
octahedron is defined in our description by the direction 
of the off-center V-ion displacement as is shown in Fig.^ 
Then, for a separate chain, four degenerate ground state 
structural configurations correspond to the following de- 
composition of the local z-axes and their variation along 
a chain: Z( = ±[bcosa± (— l)^(a x b)sina], with the 
tilt angle a ~ 45°. Because of weak interchain interac- 
tions and the influence of the side ionic groups, the exact 
degeneracy can be partly removed. We suggest, how- 
ever, that the remaining degeneracy can lead to a domain 
wall formation withing a chain, which does not contradict 
the X-ray measurements reported in P| . Chain segments 
with differently oriented bonds (domains) can sponta- 
neously form at high temperatures during the growth and 
preparation of the sample. A chain domain wall with par- 
allel neighboring local za- and zb- axes corresponds to 
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the ground-state d^y and d^y orbitals that he in paral- 
lel planes and thus have no short-parth SE connection. 
Therefore, the exchange interactions between end point 
spins of neighboring segments are strongly suppressed re- 
sulting in broken magnetic bonds which we call defects 
below. 

Chain segmentation is an attractive possibility to ex- 
plain the observed low-temperature Curie term in the 
magnetic susceptibility because segments with an odd 
number of spins lead to uncompensated S = 1/2 states. 
In the following we calculate this contribution. The ef- 
fective Hamiltonian used is the one-dimensional isotropic 
AF Heisenberg model, H = J^^S^S^+i, and we ignore 
the staggered-field terms discussed before in the following 
analysis. 

For the joint probability of having no defect between 
positions and x, but one defect at position a; on a par- 
ticular chain, or the probability to find a spin chain of 
length X, we use a generalized Poisson distribution of the 
form 

P{x)^f{x)exp(^~J^f{t)dty (4) 

Here, f{x) describes the differential probability to find 
a defect at the distance x measured from a defect lo- 
cated at the origin. For the standard Poisson distribu- 
tion, f{x) = 1/p is just a constant, i.e., is independent 
from the distance to the origin. 

A Poissonian distribution of chain lengths overempha- 
sizes the number of short chains. Instead, we expect that 
their fraction is small and the distribution peaks at a fi- 
nite value because defects which are close by "repel" each 
other due to their larger elastic energies. Lacking a de- 
tailed microscopic model for this mechanism, we simply 
assume that the probability of defects at a distance x 
grows linearly with x like f{x) = f-p^x, leading to 

= |p2^exp(-JpV) . (5) 

Here, p is the number of broken bonds per unit length of 
an infinite chain, i.e., the inverse average chain length. 
This is the Wigner distribution, well-known from quan- 
tum statistics. 



VI. NUMERICAL PROCEDURE AND RESULTS 

We use the finite-temperature Lanczos method 8] to 
calculate the eigenvalues, eigenvectors, and thermody- 
namic properties for the segmented chains. To comply 
with the effect of broken bonds, we use open bound- 
ary conditions. Then the only symmetries remaining 
to reduce the size of the problem are [H, S^] — and 
[H, Sz] = 0. While the latter can be implemented eas- 
ily using an appropriate basis, there is no simple way 
to efficiently incorporate the former. To compute finite- 
temperature expectation values, we have to store all cal- 
culated Lanczos eigenvalues and all of the corresponding 




< 

i- 0.08 



^ -I 
0.04 -I 

0.00 1 1 1 1 1 1 1 

0.0 0.5 1.0 1.5 2.0 

KbT/J 

FIG. 3; Temperature dependence of the uniform magnetic 
susceptibility of Heisenberg chains of length 23 (solid line) 
and 24 (dotted line) with open boundary conditions. 

eigenvector coefficients in terms of the basis, which is the 
limiting computational requirement here, in contrast to 
ground-state calculations. 

We have calculated the uniform magnetic susceptibility 
Xn for Heisenberg chains with length N between two and 
24 sites, defined by 

where (. . .) denotes the thermal average, Na is the Avo- 
gadro constant, /zq the magnetic permeability, g the gy- 
romagnetic ratio, fiB the Bohr magneton, and fee the 
Boltzmann constant. For the present nonmagnetic sys- 
tem (5*°*) = 0. 

To illustrate the results. Figure O shows the tempera- 
ture dependence of the magnetic susceptibility for chains 
of length 23 and 24. According to Bonner and Fisher [T^ . 
the maximum of the susceptibility for an infinite Heisen- 
berg chain is given by x™'*'^ J/(iVA^o.9^MB) ~ 0.147 and 
is reached at a temperature of fcB?inax/'/ ~ 0.641. For 
the 24-site chain, the corresponding results are 0.153 for 
the maximum susceptibility, and 0.601 for the position of 
the maximum. 

From the low-temperature upturn of xn {N odd), we 
have extracted the effective moment fj,es{N) for each 
value of N. These data are shown in Figure 0] as the 
open squares. The effective moments as a function 
of the inverse chain length scale perfectly according to 
{N) = ypv ^B- The factor \/3 originates from 
the definition of the effective moment via Curie's law, 
Xc = NAp,ofJ,lff/{3kBT). 

The effective moment as a function of the mean chain 
length N after averaging over chains of different size 
(even and odd) is shown in Figure |3| as open circles, 
representing the Wigner distribution and diamonds, rep- 
resenting the Poisson distribution. Data are shown for 
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FIG. 4: Calculated effective moments in units of the Bohr 
magneton as a function of the inverse chain length. Squares 
denote the results for chains of a single fixed length, circles for 
the Wigner average, diamonds for the Poisson average with 
mean A'^ over chain lengths between two and 24. The line 
denotes the dependence fics{N) — y^3/N /ib- 



iV < 11, because for a larger mean length, the chains of 
size greater than 24 would contribute significantly to the 
average susceptibility. 

The main effect for both averages is to reduce the mag- 
netic moment for short mean chain lengths. For the 
Wigner distribution, this is due to the linear suppres- 
sion of the weight of the susceptibilities for the short- 
est chains, whereas for the Poisson distribution, the mo- 
ment reduction is due to a comparatively large weight 
of the longer chains in the partition sum. Independent 
of how the average is performed, the effective moment is 
reduced because chains with an even site number, and 
thus no Curie term at low T, are also included in the 
average. The data suggest that for mean chain lengths 
N larger than 10, the average effective moment depends 
only weakly on the details of defect distribution. 

In Figure |31 the experimental data for the susceptibil- 
ity X of Ba2V309 are plotted as X'T versus temperature 
T. A fit to these between T = 2 K and 10 K of the form 
X — C/T+xw with C = NAfJ-afJ-lff/ (Sfce) yields an aver- 
age moment of fics ~ 0.2/xb. A mean chain length of the 
order of 75 spins would be necessary to explain the ex- 
perimentally observed Curie upturn by the segmentation 
of VO2 chains. 



Ba2V3 09. We find that within our approximate cal- 
culation the staggered field mechanism is far too small 
to explain the observed effective moment. Therefore we 

100r 




FIG. 5; Curie plot of the susceptibility of Ba2V309. The 
straight line illustrates a fit of the form x = Co/{T + Ocw) 
for temperatures T between 200 K and 400 K. In the inset, 
the same data are plotted as susceptibility times temperature 
versus T for low temperatures. 



propose a different mechanism for the Curie term which 
is based on chain segmentation due to broken magnetic 
bonds which leads to uncompensated effective moments 
on segments with odd number of spins. We find that a 
rather low concentration of ^ 10~^ broken bonds can 
explain the observed Curie term. A definite distinc- 
tion between the two models requires the investigation 
of the field-orientation dependence of the susceptibility 
in monocrystals of Ba2V30g. 
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